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We determine the general conditions for the existence of Godel, Einstein static, and de Sitter 
universes in gravity theories derived from a Lagrangian that is an arbitrary function of the scalar 
curvature and Ricci and Riemann curvature invariants. Explicit expressions for the solutions are 
found in terms of the parameters defining the Lagrangian. We also determine the conditions on the 
Lagrangian of the theory under which time-travel is allowed in the Godel universes. 
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I. INTRODUCTION 



We consider the conditions required for gravity theories that are derived from Lagrangians that are functions of the 
scalar curvature and Ricci and Riemann curvature invariants to possess solutions which are homogeneous space-times of 
[j^ ■ the Godel, Einstein static, and de Sitter forms. In the Godel case we determine the conditions for the existence or non- 
| existence of closed time-like curves in these universes. These three homogeneous space-times, and the investigations of 
^v^j ■ their stability, have played a central role in our understanding of the dynamics of general relativity and in the possible 
astrophysical consequences of general relativistic cosmologies. Higher-order modifications of general relativity are of 
importance in assessing the corrections that might be introduced to general relativity in high-curvature environments 
and can also be of use in explaining the late-time acceleration of the universe. Furthermore, investigations of these 
theories allows for an evaluation of the special nature of general relativity itself. With this in mind we have recently 
provided a detailed analysis of the cosmological consequences of gravity theories with power-law Lagrangians Q, Q| . 
Here, we extend and broaden this study to include the conditions under which the Godel, Einstein static, and de 
I ■ Sitter universes exist in a wider class of non-linear gravity theories. 
' In section [H] we present the relevant field equations for the gravitational theories that we will be considering. In 
section ITTT1 we find the conditions for the existence of Godel universes in these theories and perform an analysis of the 
conditions for the existence of closed time-like curves in these solutions. Sections II V I and IV! contain the conditions for 
q-i the existence of Einstein static and de Sitter universes, respectively, and in section IVTl we summarize our results. 
I ' 

S-H 



II. FIELD EQUATIONS 

In this paper we consider gravitation theories derived from a function of the three possible linear and quadratic 
contractions of the Riemann cur vature tensor; R, R ab R ab and R a bcdR abcd [M- The relevant weight-zero scalar density 
. . i for this general class of theories is then given by 

C G = X - 1 V9f(X,Y,Z) (1) 

where f(X, Y, Z) is an arbitrary function of X, Y and Z which are defined by X — R, Y = R a bR ah and Z = R a bcdR abcd ] 
X is an arbitrary constant which can be determined by the appropriate Newtonian low-curvature limit. The action is 
obtained, as usual, by integrating this density together with that of the matter fields over all space. The addition of 
supplementary terms to the density in order to cancel total divergences which can be transformed to integrals on 
the boundary can be problematic (see e.g. 3) and so, for simplicity, they will all be assumed to vanish. 

Taking the first variation of the action derived from , together with that for the matter fields and a cosmological 
constant, then gives 



61 = J cM(-P ab - g ab A + ^T ab )5g ab 
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where 

P ah = -\}g ab + fxR ab + 2f Y R c{a R b) c + 2f z W d < a R b \ de + f X ;cd{g ah g cd - g ac g bd ) 

+ D(f Y R ab ) + g ab (f Y R cd )., cd - 2{f Y R c{a ) b) c - KfzR d{ab)c )-,cd. (2) 

Here, A is the cosmological constant (defined independent of f(X,Y,Z)) and T ab is the energy-momentum tensor of 
the matter. The notation denotes partial differentiation of / with respect to N. A derivation of equation (J2J is 
given in the appendix. Looking for a stationary point of the action requires setting the first variation to zero, giving 
the field equations 

Pab = TjTab ~ ffafcA. (3) 

These field equations are generally of fourth-order, with the exception of the cases in which the function / is linear 
in the second derivatives of the metric, which notably includes general relativity where f = X. This property makes 
these equations particularly difficult to solve. Considerable simplification occurs if we assume that the three curvature 
scalars X, Y and Z are constant. In this case, the expression reduces to 

P ab = -\fg ab + R ab (fx - 2Rfz) - 2R acdb R cd (f Y + 4/ z ) + X -g ab {X 2 - 4Y + Z)f z (4) 



where use has been made of the identities Q 



j^abcd ^ \Z\R ad -\- R aC d 

Tjac d 1 r> ad , riabed / > , r>a r>dc 

it ■ — :r-fi. i it ttbe ~r it ,.ri 

' c 2 ' 



R ab ;ab — TI^R 



2R a cde R bcde + 2g ab Y + 2RR ab = ^g ab Z - 4R acdb R cd + 4R a c R bc + ^g ab R 2 . 

Equation jg only of second order and is therefore a significant simplification of the original system of equations. 
It is the solutions of these equations that we will now study. 

There are a number of highly symmetric space-times in which the curvature scalars X, Y and Z take constant 
values, including the Godel [5|,[1,0j Einstein static and de Sitter universes 0. We will investigate solutions of 
the field equations specified by for these three homogeneous space-times. 



III. GODEL UNIVERSES 

The Godel universe 0> 0> * s a homogeneous space-time given by the line element 

ds 2 = -{dt + C(r)d^) 2 + D 2 {r)dijj 2 



-dt 2 - 2C{r)dtdij; + G(r)dip 2 



dr 2 
-dr 2 



dz 2 
-dz 2 



(5) 



where 



C(r) 

D{r) 
G(r) 



4f2 2 fmr\ 
— Sinh {—) 



rn 

— sinh(mr) 
m 

D 2 (r) ~C 2 {r) 

4 . 2 /mr\ 

— - sinn 

m 2 \ 2 J 



1 - 



sinh 



mr\ 
~2~) 



and m and £1 are constants. The existence of Godel universes in one particular f(X,Y) theory has been previously 
studied by Accioly where a solution was found in vacuum. We extend this analysis to the more general class of 
theories above and to universes filled with matter fluids. 
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The Godel universe is of particular theoretical interest as it allows the possibility of closed time-like curves, and 
hence time-travel 0, 0] . The condition required to avoid the existence of closed time-like curves is ITU 

G(r) > for r > 



m 2 > 4fi 2 . (6) 

By investigating the existence of Godel universes in this general class of gravity theories we will also be able to 
determine those theories in which the time-travel condition © is satisfied. Consequently we will be able to determine 
those theories in which time travel is a theoretical possibility. 
It is convenient to work in the non-holonomic basis defined by 

e<°>0 = e (2) 2 = e<- 3 \ = 1 
e(°\=C(r). 

The line-element (jSJ) then becomes 

ds 2 = -(e^Y + (#«) 2 + (e^f + (e^f 

where the one-forms 9^ are given by 6^ — e^adx 11 (capital Latin letters denote tetrad indices and lower case Latin 
letters denote space-time indices). The inverses of e~ A \ can be calculated from the relations e^ A \e a , B -> = 5^ a \b) and 
the non-zero elements of the Riemann tensor in this basis are then given by 

#(Q)(1)(0)(1) = #(0)(2)(0)(2) = ^ 2 

#(i)(2)(i)(2) = 3fi 2 - m 2 . 

The perfect-fluid energy-momentum tensor is defined in the usual way with respect to the comoving 4-velocity 
U a = (1,0,0,0) and its covariant counterpart U a — (— 1, 0, — C(r), 0) such that its non-zero components in the non- 
holonomic basis are given by 

T (o)(o) = P and T {i){i) = T {2)(2) = T (3)(3) = V- 
The field equations, ©, for this space-time can then be manipulated into the form 



= (2fl 2 - m 2 )(/ x - 4[fi 2 - m 2 ]f z ) + 2{f Y + 4/ z )(8ft 4 - 5Q 2 m 2 + m 4 ) (8) 
|(p + p) - m 2 (f x - 4[Q 2 - m 2 ]f z ) + 4(/ y + 4f z )n 2 (2Q 2 - m 2 ). (9) 

Solving the field equations has now been reduced to solving these three algebraic relations for some specified /(A, Y, Z). 



A. f = f(X) 

For p + p we see from © that fx 0. From equation (JSJ it can then be seen that m 2 = 2f2 2 , as in general 
relativity. Therefore for any theory of the type / = /(A) the inequality 10 is not satisfied and closed time-like curves 
exist, when p + p ^ 0. 

It now remains to investigate the case p + p = 0. It can immediately be seen from @ that we must have fx = 
in order for a solution to exist. This sets the relation between m and f2 and automatically satisfies equation (JSJ. The 
required value of A can then be read off from equation Q. The condition fx — will now be investigated for a 
variety of specific theories. 
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1. f = f(X)=X + aX 2 

Theories of this kind have been much studied p^lll6| as they have a number of interesting properties, not least of 
which is that they display divergences which are normalisable at the one loop level P"t|. 
In these theories, the condition fx = is equivalent to 

m 2 = Q 2 + —. 

4a 

The condition under which this theory then satisfies the inequality (jSJ, and hence does not admit closed time- like 
curves, is ft 2 ^ (12a) _1 . Therefore, for any given theory of this kind, defined only by a choice of the constant a, 
there is a range of values of Q for which closed time-like curves do not exist, when a > and p + p = 0. However, 
when a < this condition is never satisfied and closed time-like curves exist for all values of f2. 



2. f = f(X)=X+^ 

Theories of this type have generated considerable interest as they introduce cosmological effects at late times, when 
R is small, which may be able to mimic the effects of dark energy on the Hubble flow [jjj, [2(| . The square in the 
factor a 2 is introduced here as these theories require a positive value for this coefficient in order for the field equations 
to have a solution. 

When fx = we have the two possible relations 

m 2 = Q 2 ±- 
2 

where a is the positive real root of a 2 . The upper branch of this solution then allows the condition © to be satisfied 
if 6f2 2 ^ a; or, for the lower branch, if 6fl 2 ^ —a. For any particular value of a > 0, the upper branch always admits 
a range of f2 for which closed time-like curves do not exist. For the lower branch, however, the inequality @) is never 
satisfied and closed time-like curves are permitted for any value of f2. 



3. f = f(X) = \X\ 



1+6 



This scale-invariant class of theories is of interest as its particularly simple form allows a number of physically 
relevant exact solutions to be found 0, 0, El ■ In order for solutions of the Godel type to exist in these theories 
we must impose upon S the constraint 5^0. 

The condition fx=0 now gives 

m 2 = n 2 . 

Evidently, the condition © is always satisfied in this case: Godel solutions always exist and closed time-like curves 
are permitted for any value of the vorticity parameter Q. 

B. f = f(X,Y,Z) 

In this general case, equations JSJl and © can be manipulated into the form 

fx m 2 m 2 )f z = 4n2(4fla _ ro2) ! X{P + P) (10) 

(2Q 2 - m 2 ) 
^ + 4/ * = - 80W-V P + P) 
when m 2 ^ 4S1 2 . The special case m 2 = 4f2 2 gives 

P + P = (12) 
/x=4f! 2 (/ y + / z ). (13) 
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First we consider p + p > 0. It is clear that to 2 = 4ft 2 is not a solution in this case. Equations (|10|) and 
show that in order to satisfy the inequality and avoid the existence of closed time-like curves, the following two 
conditions must be satisfied simultaneously 

fx + 4(to 2 - ft 2 )/ z < (14) 
Iy + ifz < 0. (15) 

For p + p < the inequalities must be reversed in these two equations. It is now clear that if it is possible to construct 
a theory which has a solution of the Godel type without closed time-like curves for a non-zero p + p of a given sign, 
then this theory cannot ensure the non-existence of closed time-like curves for the opposite sign of p + p. 

It remains to investigate the case p+p = 0. In order to have a solution in this case we require that the two equations 

fx - 4(ft 2 - m 2 )f z = (16) 

fv +4.fz = Q (17) 

are simultaneously satisfied, for m 2 ^ 4ft 2 , or 

fx=^ 2 (f Y + f z ) (18) 

for m 2 = 4ft 2 . Any theory which satisfies l|18l) . therefore, does not allow the existence of closed time-like curves when 
p + p = 0. 

We now illustrate these considerations by example. 

1. f = X + aX 2 + j3Y + jZ 

Using the identity 

(^/g(X 2 — 4Y + Z)) = pure divergence 



it is possible to rewrite / as / = X + aX 2 + (3Y where a = a — 7 and f3 = (3 + 47. This is the theory considered by 
Accioly [13 • 

For the case p + p > 0, the inequalities l|14(l and (|15f) then become 

(3 < and a > — (19) 
2\K\ 

The first of these inequalities can be satisfied trivially. The second can be only be satisfied for all £1 if \R\ has some 
non-zero minimum value. From the field equations © and 10 we can obtain 

(2ft 2 - m 2 )(l + 4a[n 2 - m 2 ]) = -2/3(8ft 4 - 5fl 2 m 2 + m 4 ) 

which gives two solutions for to as a function of ft, a and (3. Substituting either of these values of to into R = 2(ft 2 — m 2 ) 
then gives that R — > as ft 2 — * — ^j-. It can now be seen that the second equality in 1)19(1 cannot be satisfied for all ft 

if a is finite. It is, therefore, not possible to construct a theory of this type which excludes the possibility of closed 
time- like curves for all ft, when p + p 7^ 0. 

For p + p < the inequalities in H19JI must be reversed. In this case (3 > and a > allows a range of ft in which 
closed time-like curves are not permitted and (3 > and a < does not allow closed time-like curves for any values 
of ft. 

The case p + p = was studied by Accioly [l(j where it was found that the equation H18|) is satisfied if 

4ft 2 = to 2 = 5-^, 

(3d + (3) 



hence closed time-like curves do not exist in Godel universes for these theories, when p + p = 0. 
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f(X,Y,Z) 




Additional 


Closed time-like curves exist when 




conditions 


p + p > 


P + P < 


P + P = 


X + aX 2 




Q > 


/ 


/ 


//x 






Q < 


/ 


/ 


/ 


x + — 




+ve branch 


/ 


/ 


//x 






— ve branch 


/ 


/ 


/ 


\x\ 1+s 




8 > 


/ 


/ 


/ 


A -+- q:A p i -f- 'yZr 


o: 


— 7 > U, p + 47 < U 


//x 


/ 


X 




Q 


- 7 < 0, /3 + 4 7 <0 


/ 


/ 


X 




a 


- 7 > 0, /3 + 4 7 >0 


/ 


//x 


X 




a 


- 7 < 0, /3 + 4 7 >0 


/ 


X 


X 


erf + pY 




q > 0, /3 < 


/ 


/ 


X 






a > 0, /3 > 


/ 


X 


X 






a < 0, /3 < 


X 


/ 


X 






Q < 0, p > 


/ 


/ 


X 



TABLE I: A summary of the conditions under which closed time-like curves can exist in Godel universes, for various different 
gravitational theories, defined by f(X,Y,Z): / denotes their existence for all values of Q and X denotes that they are not 
allowed for any value of £1. The symbol //x means that closed time-like curves are allowed to exist for some restricted range 
of O only; the ranges are given in the main text. 

2. f = aX + PY 

This class of theories is introduced as an example which excludes the possibility of closed time-like curves for 
p + p > 0. The inequalities l|14f> and (|15|l in this case are 

a < and [3 < 

which can be trivially satisfied. For p + p = the only non-trivial solution is given by m 2 = 4f2 2 , the value of Q then 
being given in terms of a and (3 by (|18ll . 

This example shows explicitly that it is possible to construct a theory in which closed time-like curves do not occur 
in Godel universes when p + p ^ (though we do not consider it as physically viable as a is required to have the 
'wrong' sign [T5I |23| ) . 

Table [Q summarises the results found in this section. 

IV. EINSTEIN STATIC UNIVERSES 

The Einstein static universe is a homogeneous and isotropic space-time with line element 

dr 2 

ds 2 = -dt 2 + — + r 2 (d6 2 + sin 2 6W0 2 ). 

(1 — kt z ) 

Here, k parametrizes the curvature of the space-like slices orthogonal to t and the scale-factor has been rescaled to 1. 
For a universe containing pressureless dust the field equations Q can now be written as 

XP = ^fx + 16K 2 (.fY + fz) (20) 
A = i/-2 K / x -8« 2 (/y + f z ). (21) 

It can now be seen immediately that solutions exist for an Einstein static universe for any f(X, Y, Z) that is differ- 
entiable in all its arguments. The corresponding values of p and A are simply read off from equations (|20|l and l|21|l . 
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It remains to be studied under what circumstances these solutions are stable. The investigation by Barrow, Ellis, 
Maartens, and Tsagas @ shows that this is an issue that depends upon the material content and the equation of state 
of matter in a delicate fashion. In general relativity, there is first-order stability against density perturbations when 
the sound speed exceeds a critical value (l/y/E of the speed of light) because the Jeans length exceeds the size of 
the universe |8L l24j . However, there is instability against homogeneous gravitational- wave modes of Bianchi IX type 
p|. In general, we expect that a universe with compact space sections and Killing vectors to display linearisation 
instability ,25| . The function space of general solutions to Einstein's equations possesses a conical structure at these 
particular special solutions and there are an infinite number of perturbation expansions tangential to the conical point 
which do not converge to true solutions to the field equations. This is only ensured if further constraints are satisfied 
and some investigation of these problems was made by Losic and Unruh |2q . Similar studies could be performed for 
the new family of solutions that we have identified here. 



V. DE SITTER UNIVERSES 



The line element for the maximally symmetric de Sitter vacuum universe can be written as 

ds 2 = -dt 2 + e 2 ^\dr 2 + r 2 (dQ 2 + sin 2 6d<j) 2 )) 

where A is the cosmological constant. For de Sitter space-time all components of the Riemann and Ricci tensors can 
be written in terms of the Ricci scalar using the equations 

1 

Rabcd — —-^R{gad9bc — 5oc56d) 
1 

Rab — -^Rgab 

where the Ricci scalar is R = 4A. The field equations J3J) can now be reduced to the single equation 

\f-A = Af x + 2A 2 / y + ^A 2 f z . (22) 

This equation must be satisfied by / if the de Sitter universe is to be a solution in any particular gravitational theory. 
This problem reduces to that studied by Barrow and Ottewill [15| in the case / = f(X). This result establishes the 
situations where inflation of de Sitter sort can arise from higher-order corrections to the gravitational lagrangian. 
In the case where / = f(X ) al one it is appreciated that the resulting theory is conformally equivalent to general 
relativity plus a scalar field 27, 28] with an asymmetric exponential potential and so either de Sitter or power-law 
inflation is possible. Our results establish when de Sitter inflation is possible in situations where the other invariants, 
Y and Z contribute to the Lagrangian, and the conformal equivalence with general relativity is broken. The stability 
of the de Sitter solutions in these theories will be studied elsewhere. 



VI. CONCLUSIONS 



We have analysed the existence conditions for solutions of Godel, Einstein, and de Sitter type to exist in gravity 
theories derived from a Lagrangian that is an arbitrary function of the curvature invariants R, R a bR ab and R a i, c dR abcd . 
The existence conditions have been systematically explored and a number of special choices of the Lagrangian function 
that are of physical interest were worked out explicitly to display the form of the solutions in terms of the Lagrangian 
properties. In the Godel case there is a simple condition on the metric parameters which reveals whether or not time- 
travelling paths exist in the space-time. We have evaluated this condition in general and explicitly in the examples 
studied. A range of situations exist in which time-travel is either possible (as in general relativity) or impossible 
depending on the form of the gravitational Lagrangian function defining the theory. We find that it is not possible 
to construct a theory, within this class, for which the existence of closed time-like curves is forbidden for all perfect 
fluids. The conditions for Einstein static and de Sitter solutions were also found. 
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VIII. APPENDIX: DERIVATION OF EQUATION © 

The variation of the action derived from integrating the density JQ) over all space is 
SI G = X' 1 J dn^g[^fg ab Sg ab + f x SX + f Y SY + f z SZ 
dn^\\fg ab 5g ab - fx(R ab S 9 ab - g ab SR ab ) 



where use has been made of 



and 



2f Y (R< a R b \Sg ab - R ab SR ab ) ~ 2f z (R cde (b R^ edc Sg ab - R a bcd 6 R a bcd )] (23) 



5g ab = -g ac g bd 5g cd 



s Vg = ^g ab Sgab- 



The subscript a here denotes that we are considering the gravitational part of the action only. Using the relations Q 

1 



SR ab - ~-g c (Sg ab;cd + Sg cd , ab - 6g ac , bd - 5g bd;ac ) 



and 



we can then write 



$R a bcd = -^9 ae {5g ed - bc + Sg eb - dc - Sg db - 



ec u yec 



bd — Sgeb-.cd + Sg cb - ed ) 



fxg^SRab =S ~fx;cd(9 ab 9 Cd - 9 aC 9 bd )S 9 ab 



2f Y R ab SR a , 



n(f Y R ab ) + {f Y R cd ), cd g ab - 2(f Y R< a ), \ 5g ab 



2fzR a bcd SR a bcd * 4(fzR c{ab)d );cd6gab 

where ~ means equal up to terms which are pure divergences. Such terms are irrelevant here as they can be transformed 
via Gauss's theorem to terms on the boundary which are assumed to vanish. Substituting these expressions back into 
(|23|l and making the definition 

SI G = -X' 1 I dnP ab 6g ab 



then gives equation |2"|. which completes the derivation. 
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